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1 . INTRODUCTION 


1.1 OPS BENCH PROGRAM DESCRIPTION 

The OPS Bench Navigation Program is a single-vehicle weighted least squares 
orbit detev’mination (OD) program with the added capability of generating and 
storing ephemerides for three Tracking Data Relay Satellites (TORS). This 
program was developed under NASA Contract NAS9-15722 for use in verification 
of the OPS MCC Ground Navigation Program. 

Both ground-based and satellite-relayed observations may be processed by 
the OPS Bench Program. The ground-based observations may be in-plane and out- 
of-plane angles, range, or Doppler. The satellite-relayed observations may be 
range or Doppler. 

The ICOASV predictor in the OPS Bench Program is designed to propagate a 
vehicle state vector in free flight under the influence of significant 
natural forces in the solar system. The scope of the modeling is limited to 
perturbative forces consistent with free-flight techniques. On option, the 
models for atmospheric drag, solar radiation, and venting may be included in 
the predictor computations. Integration is performed in double precision using 
a variable order Gauss-Jackson numerical integration scheme. The normal inte- 
gration mode is the Encke-beta mode, but it is possible to select via user 
input either the Encke-time mode or the Cowell mode. 

The OPS Bench Program employs a differential correction processor to solve 
the weighted least squares problem and update the vehicle state and error 
covariance matrix. Parameters which may be included in the differential correct! 
process are components of vent forces, drag multiplier, and relay biases. 

Output of the OPS Bench Program in the execution of an OD solution consist 
of a listing of input data, including a full listing of data batches employed 


in the OD process, a summary of vehicle trajectory, and a complete summary of 
the differential correction iterations. Vehicle trajectory, when such a 
trajectory is required, may be displayed in a number of coordinate frames. 

The OPS Bench Program is written in double-precision FORTRAN for the 
UNIVAC 1108 - EXEC 8 operating system. 


1.2 USER'S MANUAL SUMMARY 


The OPS Bench Program User's Manual describes in detail all information 
required by the user to operate successfully the OPS Bench Program. 

Section 2 discusses the tape and card inputs to the OPS Bench Program, and 
complete details of the FORTRAN NAMELIST input format are provided. 

Section 3 offers complete descriptions of each of the namelist entities, 
and a listing of inputs available to the user in each namelist. 

Section 4 provides input tape unit assignments and program deck set up for 
a Bench Program run. 

Section 5 contains sample cases which illustrate the output of typical OPS 


Bench Program runs. 


2. PROGRAM INPUT 


In this sectioHj the various tape and card inputs to the OPS Bench Program 
are discussed. For reference, a simplified input/output scheme for the orbit 
determination process of the Bench Program is shown in Figure 2-1 on the 







2.1 TAPE INPUT 


Tape inputs to the OPS Bench Program fall into several categories. These 
include a program tape, observation data batch tape, station characteristics 
tape, and planetary ephemeris tape. The basic format of each of these tapes 
is described in a separate subsection below. 

2.1.1 OPS Bench Program Tape 

The OPS Bench Program Tape is required for the execution of an OD solution. 
This is a baseline tape which contains the actual Bench Program software, and 
several programmer reference files. The files contained in the tape are as 
follows: 

(a) MAINTAIN File . This file consists of canned (3ADD elements wh’.h perform 
a number of functions essential to the creation and use of the Baseline 
tape. 

(b) SOURCE File . This file contains the FORTRAN source for the Bench Pro- 
gram and its utility processors. In addition, SOURCE contains the 
COMGEN Data Base. 

(c) RELOC File . This file contains the relocatable elements created by 
'the compilation of the SOURCE file elements. 

(d) ABS File . This file contains the absolute executable elements for the 
Bench Program and its utility processors. 

(e) MAPS File . The MAPS file contains the map directives necessary to 
create the ABS file elements from the relocatable elements in RELOC. 

(f) TCASES File . This fi le provides the programmer with input data for 
a set of standard test cases developed for the Bench Program. 
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1 (g) ALTERS File . This file contains the modifications used to create 

the current Baseline tape from its predecessor. 

2.1.2 Observation Data Batch Tape 

The Observation Data Batch Tape is required input for an OD in which 

{ 

j ground-tracking or TDRS-relayed data is to be processed. 

The Data Batch Tape consists of an arbitrary number of batches of navi- 
‘ gation tracking data, each batch accompanied by a batch header array of infor- 

. mation. 

: '■ The data batches occur in chronological order based upon the time of the 

■I first valid data observation contained in the batch. Each data batch consists 

of a time sequence of data frames, each containing the following information: 

” (a) Data frame time tag 

' . (b) Angle observables (XY30,XY85,Az,El) 

(c) Range observable 

[ (d) Doppler observable 

: (e) Doppler count interval 

L 

b , Each data batch contains a maximum of 100 data frames. 

, Each, data batch is preceded by its own batch header. This header is an 

; array containing the following information: 

i (a) NR = Receiver ID 

NX = Transmitter ID 

: (b) NTl = TDRS Forward Link ID 

NT2 = TDRS Return Link ID 
^ (c) NV = Vehicle ID 

DT = Data Tape 
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(d) NF = Number of Data Frames 

(e) = Transmitter Direct Frequency 

(f) = Receiver Direct (Reference) Frequency 

(g) “ Return TDRS Translation Frequency 

r 

(h) = Doppler Model Frequency Multiplier 

(i) = Doppler Bias Frequency 

(j) A|^ = Range Ambiguity Interval 

I'.* ’ • 

(k) BT = Batch Time 

■ t 

2.1.3 Station Characteristics Tape 

The Station Characteristics Tape is required input for an OD problem. 

This tape provides the identifying characteristics of various STDN tracking 
stations currently employed for ground navigation of the Shuttle Orbiter. The 
information supplied for each tracking station is the following: 

STA(l) = Station ID (JSC acronym) 

STA(2) = Geodetic Longitude 
STA(3) = Geodetic Latitude (GDLAT) 

STA(4) = Geocentric Latitude (GCLAT) 

STA(5) = sin (GDLAT) 

STA(6) = cos (GDLAT) 

STA(7) = sin (GCLAT) 

STA(8) = cos (GCLAT) 

t 

STA(9) = Altitude above Ellipsoid 
STA(IO) = Station Radius (R^) 

STA(ll) - STA(22) = Refraction Moduli (N^) 


STA(23) - STA(34) = Atmospheric Scale Heights (Hg) 

STA(35) = Ellipsoid Radius (R^) 

STA(36) = Rg[STA(8)] + [STA(9) ]• [STA(6) ] 

STA(37) = Rg[STA(7)] + [STA(9)]*[STA(5)] 

STA(38) = 2-way Doppler Coefficient (Wq 2 ) 

STA(39) = 3-way Doppler Coefficient 
STA(40) = Range Coefficient (Wj^) 

STA(41) = Angle Coefficient (W^) 

2.1.4 Planetary Ephemeris Tape 

The Planetary Ephemeris Tape is required for the calculation of vehicle 
acceleration due to the planetary disturbing bodies selected in the gravita- 
tional potential model of the ICOASV predictor. The current tape in use is 
the PE50D tape generated from the JPL DE19 series ephemeris tapes. It provides 
the ephemerides of all planets and the Moon, spanning the time December 1961 to 
December 1994. 

The time interval for planetary ephemeris points is 4 days, with the 
exception of Mercury, which has a 2 day interval. The interval between the 
Moon's ephemeris points is 12 hours. The planet data consist of heliocentric 
position and velocity components with their modified second and fourth differ- 
ences. Data for the Moon consist of geocentric position and velocity components 
and the associated modified second and fourth differences. Units for all 
position and veloci ty- components are ER and ER/hr. All data are referenced 
to the M50 coordinate system and are represented in double precision. 


2.2 CARD INPUT 


The FORTRAN NAMELIST input scheme is employed by the OPS Bench Program 
for card input data. The NAMELIST scheme prc/iries a streamlined method by 
which data can be input according to certain oataijories called namelists. In 
each namelist, inputs are identified by variable names, and data is input in 
a free-field format. 

2.2.1 Bench Program Namelists 

Input to the OPS Bench Program by cards is divided into thirteen separate 
namelists. Each namelist provides a set of data variables which collectively 
furnish data for a specific function. For example, the namelist MISVEH allows 
the user to input the data necessary to describe the particular Mission and 
VEHicle characteristics for the 00. 

The following is a complete list of the namelists in the Bench Program 
input scheme: 

(s) SYSTEM - System Parameters 

(b) MISVEH - Mission and Vehicle Characteristics 

(c) VADRAG - Variable Drag Model Parameters 


(d)' 

PRBCTL 

- OD Problem Control 

(e) 

SBCTL 

- Superbatch Control 

(f) 

BTBCTL 

“ Batch-to-Batch Control 

(g) 

TRAJ 

- Trajectory Control 

(h) 

ATOEDT 

- Automatic Edit Parameters 

(i) 

SOLFOR 

- Solve-for Parameters 

(j) 

COVARS 

- A Priori Covariances 

(k) 

MANEDT 

- Manual Editing Options 


(l) BRESPL - Batch Residual Plotter Control 

(m) BIASCT - Observation Data Bias Control 

2.2.2 Namelist Format 

The general format for the Bench Program namelist input for a single OD 

case is as follows: 

Column 2 
I 

$ NAMEl 

! (NAMEl Data Cards) 

$ END 
$ NAME2 

(NAME2 Data Cards) 

$ END 
$ NAMEj 

(NAMEj Data Cards) 

$ END 
$ ENDCAS 

A namelist that is called for data input is opened with a card $NAME having a 
$ in Column 2 followed by the formal name for that namelist. Following the 


namelist card are the input data cards for that namelist. These data cards 
will follow format specifications detailed in Section 2.2.3 below- The name- 
list is then closed with the $END card; as illustrated, the $ appears in 
Column 2, followed by the word END. The cards for each namelist requiring 
input follow this format. These namelist cards are stacked as shown above into 
a single namelist input deck which is closed with the $ENDCAS card. The order 
of appearance of the namelists is unimportant, and is left to the user's dis- 
cretion. The $ENDCAS card terminates the user input for the particular Bench 
Program case. 

If more than one Bench Program case is to be submitted in a single 
run, the namelist input decks for each problem appear consecutively in the 
program run deck, the $ENDCAS cards separating the successive namelist input 
decks. Note : The Program does not carry over input from one case to the 

next. The namelist input deck provided for each case must therefore be self- 
contained. 

One final item concerning the Bench Program's namelist format must be 
clarified. This item concerns the requirement for calling a namelist in the 
namelist input deck. All data variables available in the Bench Program are 
assigned default values. These default values may either be zeros or nonzero 
operating values. The user may elect to accept the default value for any data 
variable by simply omitting to input any data for that variable. Hence, the 
list of data variables called for input in a namelist may either be a full or 
partial list of all the data variables available in that namelist. If the full 
set of default values for a particular namelist is acceptable to the user, or 
no input for that namelist is required for the case submitted, the user is 
not required to call that namelist. For example, the SYSTEM namelist is 


provided a complete set of default values for system parameters and integrator 
options which the user may opt for by omitting to call SYSTEM in the namelist 
input deck. An example of a namelist which may not require input for an OD 
is the namelist BIASCT : if no data biases are to be applied to the observation 
data batches to be processed in the current OD, then the user omits calling 
of BIASCT . 

All default values for the data variables available in the Bench Program 
are listed in the user input descriptions of Section 3. 

2.2.3 Input Variable Format 

At the end of this subsection, an alphabetic listing of all data variables 
available in the OPS Bench Program is provided in Table 2.1. Corresponding 
to each data variable in the table there is given the namelist, or namelists, 

'to which the data variable belongs. The user may then refer to Section 3 for 
complete description of data variables. 

Input values for any data variable may take any of the following forms: 

• Integer constants, 

• Real Constants. These are written with a decimal point, and, optionally, 
they are written with an exponent consisting of E, +, or ~, or a 
combination of E and sign. 

• Double precision constants. These are written with a decimal point 
and nine or more digits, and, optionally, they are written with an 
exponent consisting of D, +, or or a combination of D and sign. 

If the input value contains fewer than nine digits, the D form must 
be used with or without an exponent (e.g., l.ODO). 

• Logical constants. These are written as T, .TRUE., F, or .FALSE.. 


• Hollerith constants. These are written as nHhhh...h, where hhh...h is 
a string of n alphanumeric characters including blanks. Six characters 
can be stored in one location. If less than six characters remain, they 
are stored left-justified with the rest of the computer word filled 
out with blanks. 

To assign input to a particular data variable, the namelist to which that 
variable corresponds must first be called, as detailed in the previous section. 
Input data cards follow this format line: Variable name can begin on a card 

anywhere execept Column 1. Since the $ cards begin in Column 2, it is suggested 
that variable names be placed beginning in Column 3 for ease in identifying 
separate namelists in the printout of user input data. The variable name is 
followed by an equal sign, which is followed in turn by the input data. In the 
case of input to a variable array, input elements are separated by commas. 

Arrays of input elements may overflow to one or more other cards, with the 
provision that no blanks should be embedded in an input constant. Blanks may be 
freely used elsewhere. 

The forms that the input data may take are the following: 

(a) VARIABLE NAME = constant . In this case, the variable called is a 
simple variable. Examples of this type of input are 

NTDRS = 2; and 

LEDIT = .TRUE. 

The first example is that of the namelist MISVEH input which specifies 
the number of TORS vehicles involved in the OD problem. The second 
example is that of the namelist BTBCTL input which specifies that the 
automatic edit option in the batch-to-batch processing mode is desired. 


(b) SUBSCRIPTED VARIABLE = constant . In this case the variable called is 


a single element of a variable array. An example of this type of 
input is 

POTENT(ll) = -2.15D-6 

This example is one in which the coefficient J 3 of the Earth geo- 
potential model, located in the eleventh position in the 600-element 
array POTENT, is given the value -2.15D-6 by the user, rather than 
allowing use of the default value -2.56D-6. POTENT is a SY STEM name- 
list input, so the usual namelist call is required for this input, 

(c) SUBSCRIPTED VARIABLE = set of constants . In this case, an ordered 
string of elements in a variable array is called for user input. 

The variable array name is followed immediately on the data card by 
the subscript of the array element at which the input is to begin. 
Following the equal sign, the constants appear in the order in which 
input is desired, and the number of constants input must not exceed 
the space available in the array following the starting element. 

Input may be individual constants, or a repeated constant. The latter 
will be of the form R*constant where R is an unsigned integer. 

Examples of this type of input are 

POTENT(IO) = 1079. 2D-6, -2.15D-6, -1.6D-6; and 

P0TENT(296) = 21*0. DO 

In the first example, the coefficients of the Earth geopotential model 
^^4 have been assigned values by the user. occupies the 
array element POTENT(IO). In the second example, the Moon geopotential 
' coefficients through C^ 3 have been assigned zeros. In actuality, 




J! ■ ■ 


only the coefficients J 2 , J 3 , C 2 £’*"3 1 ’ ^3 3 P0TENT(297), 

P0TENT(298), P0TENT(313), P0TENT(314), and POTENT(316)) have nonzero 
default values, but the method of giving them zero values shown in 
the example is easier for the user than handling each coefficient 
as a separate subscripted variable. 

In the case of a doubly subscripted variable array, an ordered 
string of elements is called for user input by columns . An example 
of this type of usage is the following: 

VTF0RC(1,1) = 20.0, -2.6, 0.8 

The VTFORC array in the SOLFOR namelist is used to input the a priori 
force levels for a maximum of three solve-for vents. Each a priori 
vent is resolved into its component forces, in units of lbs, along 
each of the vehicle body coordinate axes. For the J-th such vent, 
the XgY"> '''by"’ forces are loaded into the J-th 

column of VTFORC (into the (1,J), (2,J), and (3,J) entries, respec- 
tively). In the example above, the user has input the first a priori 
vent by components into the first column of VTFORC. Input of a second 
a priori vent could be accomplished in either of two ways, at the user's 
conveni ence: 

(1) VTF0RC(1,1) = 20.0, -2.6, 0.8 
VTF0RC(1,2) = 2.0, 0.0, 1.8 
or 

(2) VTF0RC(1,1) = 20.0, -2.6, 0.8, 2.0, 0.0, 1.8 

(d) ARRAY NAME = set of constants . In this case, the variable called, is 
an entire array. The set of input values may be individual constants, 
or a repeated constant. The set of input values specified will be 
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loaded by columns into the array beginning at the first position: 
NAME(l) for a singly subscripted array, or NAME(1,1) for a doubly 
subscripted array. The number of input values must not exceed the 
storage available in the array. 

An example of this type of array input is in input of a vehicle 
state vector in M50 coordinates (x,y,z,x,y,z) via the INVEC variable 
array in the MISVEH namelist. For convenience, this might be input 
on two cards: 

INVEC = .828614678, .311808358, .532416856, 

-.596813269, 4.097290949, -1.459561936 


Table 2.1. OPS Bench Program Input Variables 


fe, ' ■ 


VARIABLE 


NAMELIST 


** 3 ** 


BCONV 

BIAS 

BIASIN 

BIASV 

BTAB 


CDA 

CDETUT 

CDF 

CDN 

CDRAG 

CDS 

CPN 

CSA 


DA 

DAYS 

DCONV 

DELTAT 

DELTAl 

DELTA2 

DKFACT 

DKl 

DK2 


EDLIM 

EPSQR 

ETAB 


FKGAMA 


GET 


SBCTL 

SOLFOR 

BIASCT 

BIASCT 

MISVEH 


★*Q** 


VADRAG 

MISVEH 

VADRAG 

VADRAG 

MISVEH 

VADRAG 

VADRAG 

MISVEH 


*★[)** 


SBCTL 
MISVEH 
SBCTL 
TRAJ 
ATOEDT 
ATOEDT 
SOLFOR 
ATOEDT 
1 ATOEDT 


**£** 


HMS 


ATOEDT 

SYSTEM 

MISVEH 


I COVARS 

**Q** 

I MISVEH 

**l-l** 

I MISVEH 


VARIABLE 


NAMELIST 


IBATT 

ICOV 

ID 

IDCM 
IDIS 
I EXBAT 
IG 

I MODE 

INVEC 

IPRINT 

IRKST 

ISLVBI 

ISLVDR 

ITRF 

lUNITS 

IVARFX 


JDRAG 

JSOLAR 

JVENT 


KEYS 

KPR 

KREF 


LANCHL 

LAPCOV 

LEDIT 

LEXCCA 

LEXCCE 

LEXCCR 

LEXCHD 

LEXCHR 

LEXCLC 

LEXCLR 

LEXCLS 

LEXCRD 

LEXCRR 




MISVEH 

COVARS 

ATOEDT 

PRBCTL 

COVARS 

MANEDT 

ATOEDT 

SYSTEM 

MISVEH 

SYSTEM 

SYSTEM 

SBCTL 

SBCTL 

MISVEH 

BRESPL 

SYSTEM 


**j*:fc 


PRBCTL 

PRBCTL 

PRBCTL 




MISVEH 

PRBCTL 

PRBCTL 


**!_** 


SBCTL 

BTBCTL, 

SBCTL 

BTBCTL 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 
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Program Input Variables (Concluded) 


VARIABLE 


namelist 


VARIABLE 


NAMELIST 


LEXCSX 

LEXCSY 

LEXCSR 

LEXCS2 

LEXCS3 

LFORCE 

LTPROP 


MAXDV 

MAX IT 
MEDIT 
MINCA 

MINEL 

MINI 

MIN2 

MONTH 

MXELOP 


NFBAT 

NLBAT 

NORD 

NPBI 

NTDRS 

NVNTSL 


POTENT 

PRODES 


RCONV 


MAN EOT 

MANEDT 

MANEDT 

MANEDT 

MANEDT 

BTBCTL, 

SBCTL 

SBCTL 


**(v|** 


BTBCTL, 

SBCTL 

BTBCTL, 

SBCTL 

MANEDT 

BTBCTL, 

SBCTL 

BTBCTL, 

SBCTL 

ATOEDT 

ATOEDT 

MISVEH 

ATOEDT 


**[\|** 


BTBCTL, 

SBCTL 

BTBCTL, 

SBCTL 

SYSTEM 

COVARS 

MISVEH 

BTBCTL, 

SBCTL 


**p** 


SYSTEM 

PRBCTL 


**p** 


BTBCTL, 

SBCTL 


TDRS 

TEDELT 

TSTART 

TSTOP 


VCONV 

VEDELT 

VENTIN 

VTCONV 

VTFORC 

VTSTRT 

VTSTOP 


WT 

WTAB 


YEAR 


RMAX 

ATOEDT 

RMIN 

ATOEDT 

RMNEL 

(BTBCTL, 
1 SBCTL 

RPLOT 

BRESPL 

**5'A-* 

SCVBH 

COVARS 

SCVB23 

COVARS 

SCVNT 

COVARS 

SDRAG 

COVARS 

SRCHl 

ATOEDT 

SRCH2 

ATOEDT 

SSCOV 

COVARS 




MISVEH 

SYSTEM 

TRAJ 

TRAJ 




BTBCTL, 

SBCTL 

SYSTEM 

MISVEH 

BTBCTL, 

SBCTL 

SOLFOR 

S0LFOR 

SOLFOR 


**[ij** 


MISVEH 

MISVEH 


**Y** 

I MISVEH 
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3. USER INPUT DESCRIPTION 


This section provides a delineation of user inputs for the OPS Bench 
Program. Each namelist is detailed in a separate subsection: the namelist 

itself is described, and the input variables available in that namelist are 
catalogued. Included at the end of each namelist subsection is an alphabetic 
listing of the namelist variables with full description of each variable 
(name, dimension, type, capsule description, and default value). 

Provided in Figure 3.1 on the following page for the user's convenience 
is a "walk-through" flow diagram of the namelist input scheme. Each namelist 
is labeled either as optional or mandatory. If a namelist is mandatory., the 
OD requires input in that namelist. If a namelist is optional, the user may 
provide inputs or else accept the default values assigned to the variables 
in that namelist. This diagram will enable the user to ensure that sufficient 
data has been provided for the OD. 
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Figure 3.1. OPS Bench Program Input Flow Diagram (Continued) 













Figure 3.1. OPS Bench Program Input Flow Diagram (Concluded) 






3.1 SYSTEM NAMELIST 


Through the SYSTEM namelist of inputs, the user exercises control of 
parameters governing the performance of the ICOASV free-flight predictor. 
SYSTEM inputs include major integrator options, vehicle ephemeris time param- 
eters, and coefficients for the disturbing body gravitational potential 
models. 

3.1.1 Integrator Options 

Numerical integration is performed in double precision using a Gauss- 
Jackson backward difference algorithm with a Runge-Kutta starter. Integrator 
options available to the user are the following: 

(a) Integrator mode (IMODE) 

(b) Rectification tolerance (EPSQR: Encke modes only) 

(c) Integration step control (IVARFX) 

(d) Order of Gauss-Jackson backward difference (NORD) 

(e) Number of Runge-Kutta steps per Gauss-Jackson (IRKST) 

(f) Integrator printout control (IPRINT) 

(g) Shuttle ephemeris storage time step (VEDELT) 

(h) TORS ephemeris storage time step (TEDELT) 

3.1.2 Gravitational Potential Model Parameters 

ICOASV performs computation of vehicle acceleration due to the gravita- 
tional attraction of selected disturbing bodies in the solar system. (See 
Section 3,1.3, Note 1). The coefficients of the gravitational potential model 
for each of the selected bodies are input via the POTENT array. 


For the OPS Bench Program, Earth is selected as the primary body, and 
the Sun and Moon are selected as secondary disturbing bodies. Other planets 
and natural satellites are available for inclusion as disturbing bodies, but 
require program modification. (See Section 3.1.3, Note 2). 

POTENT is a 600-element array. Coefficients for the potential models of 
certain planets and the Moon can be input as follows: 

(a) Coefficients for Venus model (P0TENT(l-8)) 

(b) Coefficients for Earth model (P0TENT(9-263)) 

(c) Coefficients for Mars model (P0TENT(264-271)) 

(d) Coefficients for Jupiter model (P0TENT(272-279) ) 

(e) Coefficients for Saturn model (POTENT (280-287)) 

(f) Coefficients for Neptune model (P0TENT(288-295) ) 

(g) Coefficients for Moon model (P0TENT(296-550) ) 

(h) P0TENT(550-600) (not used) 

3.1.3 Notes to User 

Note 1: Four sources of acceleration on a vehicle in free flight can be 

modeled in the Bench Program: gravitation, atmospheric drag, solar radiation, 

and venting. Parameters for drag and vent models can be found in the name- 
list VADRAG and MISVEH . The switching flags for the drag, solar radiation,, 
and vent models are set via the problem control namelist PRBCTL . 

Note 2: The inclusion flags for planets and certain natural satellites 

are located in the IPLNET and ISATEL arrays in the 6LKDAT subroutine. 
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Tables.!. SYSTEM Namelist 






S 


t -• 


VARIABLE 

DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

EPSQR 

1 

DP 

Rectification tolerance for Encke integration 
modes 


l.OD-4 

1 I MODE 

i 

i 

1 

I 

Integration mode flag 
=0, Encke-beta mode 
=1, Encke- time mode 
=2, Cowell mode 


0 

! IPRINT 

i 

ro 

1 

I 

Integrator printout flag 
=-l, print initial and final vehicle state 
=0, no state vector print 
>0, print every N-th point computed along the 
trajectory (where N = IPRINT) including 
initial and final vehicle state 


-1 

i IRKST 

1 

I 

Number of Runge-Kutta steps per Gauss-Jaekson 
step 


1 

IVARFX 

1 

1 

I 

Indicator flag for variable or fixed inte- 
grator step size 
=0, fixed step size 
=1, variable step size 


0 

NORD 

1 

I 

Highest order backward difference in the Gauss- 
Jackson algorithm (NORD £ 15) 


8 

POTENT 

i 

600 

DP 

Coefficients for the non-spherical gravitation 
potential models of the planets and natural 
satellites 



\ 

1 

> 

•1 

i 

-i 




\ 



t,x --n: 





2*’’***l?!5! 






} 

I 

t 

i 

i 

■ f 

i 


1 

I 

i 

I 


[i 

f; 


;ri 

H 


i'i 

i\ 


Table 3.1. SYSTEM Namelist (Continued) 
VARIABLE dim TYPE DESCRIPTION 


UNITS DEFAULT 


POTENT(I) (1=1,2 8 ) 

Coefficients for Venus Model 


POTENT(l-2) 

J-j through J 2 

2*0. DO 

POTENT(3) 

^ 1 ,1 

O.DO 

POTENT (4- 5) 

C^ ^ through C^ 2 

2*0. DO 

P0TENT(6) 

^ 1,1 

O.DO 

P0TENT(7-8) 

S 2 J through 2 

2*0. DO 

POTENT(I) (1=9,10 263) 

Coefficients for Earth Model 


P0TENT(9) 

'^l 

O.DO 

POTENT(IO) 

^2 

1082.7D-6 

POTENT(ll) 

^3 

-2.56D-6 

P0TENT(12) 

J4 

-1.58D-6 

POTENT( 13-23) 

Jg through J-jg 

11*0. DO 

POTENT (24) 

'1,1 

O.DO 

P0TENT(25) 

^ 2,1 

O.DO 

POTENT (26) 

*^ 2,2 

1.57D-6 

POTENT (27-29) 

C^ 1 through C^ 3 

3*0. DO 

POTENT (30- 33) 

^ 4^7 through C^’^ 

4*0. DO 

POTENT( 34-38) 

Cg ^ through Cg g 

5*0. DO 

POTENT( 39-44) 

Ce.’i 

6*0. DO 

POTENT (45- 51) 

Cj ^ through Cj ^ 

7*0. DO 

POTENT (52-59) 

Cg’, through Cg^g 

8*0. DO 

POTENT( 60-68) 

Cg^^ through Cg^g 

9*0. DO 


# 

1 . 

I 

( 

f 


I 


Table 3.1. SYSTEM Namelist (Continued) 





ro 

00 


VARIABLE 

DIM TYPE DESCRIPTION 

UNITS DEFAULT 

P0TENT(69-78) 

Cioj through 

10*0. DO 

POTENT( 79-89) 

^11,1 ^11,11 

11*0. DO 

P0TENT(90-100) 

Cl 2 1 through C-j^ ^2 

12*0. DO 

POTENT (102- 11 4) 

^13,1 C-|3j3 

13*0. DO 

P0TENT(115-128) 

^14,1 

14*0. DO 

POTENT( 129-1 43) 

^15,1 C-J5J5 

15*0. DO 

P0TENT(144) 

^ 1 ,1 

O.DO 

POTENT (145) 

^ 2,1 

O.DO 

POTENT (146) 

^ 2,2 

-.897D-6 

POTENT (147- 149) 

S 3 1 through S 3 3 

3*0. DO 

POTENT( 150-1 53) 

S^^^ through S^^^ 

4*0. DO 

POTENT( 154-1 58) 

Sg’, through Sg’j 

5*0. DO 

POTENT( 159-1 64) 

Sg_, through Sg _5 

6*0. DO 

POTENT (165- 171) 

Sj 1 through Sj -j 

7*0. DO 

P0TENT(172-179) 

Sg’, through Sg’g 

8*0. DO 

POTENT (180- 188) 

Sg^^ through Sg^g 

9*0. DO 

POTENT( 189-198) 

$ 10,1 ’^*''■'’''9'' $ 10.10 

10*0. DO 

POTENT( 199-209) 

^ 11,1 

11*0. DO 

P0TENT(21 0-221) 

S-jp ^ through S ^|2 ^2 

12*0. DO 

POTENT( 222-234) 

$13,1 $13,13 

13*0. DO 

P0TENT(235-248) 

$14,1 ‘'"■“'^Sh 

14*0. DO 

P0TENT(249-263) 

^15,1 t^'"°^'9h S-,5j5 

15*0. DO 


' ;fy-- - - . viWi l iV - y— 




Table 3.1. SYSTEM Namelist (Continued) 


VARIABLE 


DIM type 


DESCRIPTION 


ro 

i.£) 


POTENT(I) (1=264,265 271) 

POTENT (264) 

POTENT (265) 

P0TENT(266) 

P0TENT( 267-268) 

POTENT (269) 

POTENT (270-271) 

POTENT(I) (1=272,273,... ,279) 

POTENT (272) 

P0TENT(273) 

P0TENT(274) 

P0TENT(275-276) 

POTENT (277) 

POTENT (278-279) 

POTENT(I) (1=280, 281,..., 287) 

POTENT (280) 

P0TENT(281) 

POTENT (282) 

P0TENT( 283-284) 

POTENT (285) 


Coefficients for Mars Model 

^1 

J 2 

^ 1,1 

C 2 1 through C 2 2 

^ 1,1 

$2 1 through $2 2 
Coefficient for Jupiter Model 

^2 

^ 1,1 

C 2 1 through C 2 2 

^ 1,1 

$2 1 through $2 2 
Coefficients for Saturn Model 


1,1 

C 2 ■] through C 2 2 

^ 1,1 


UNITS DEFAULT 


O.DO 
.197D-2 
O.DO 
2*0. DO 
O.DO 
2*0. DO 


O.DO 
-.296D-1 
O.DO 
2*0. DO 
O.DO 
2*0. DO 


O.DO 
.27D-1 
O.DO 
2*0. DO 
O.DO 




Table 3.1. SYSTEM Namelist (Continued) 


VARIABLE DIM TYPE DESCRIPTION UNITS DEFAULT 


POTENT (286- 287) 

S^ 1 through S^ 2 

2*0. DO 

POTENT(I) (1=288, 289,..., 295) 

Coefficients for Neptune Model 


POTENT (288) 

>^1 

O.DO 

POTENT (289) 

J 2 

.35D-2 

POTENT (290) 

^1,1 

O.DO 

P0TENT(291-292) 

C 2 1 through C 2 2 

2*0. DO 

P0TENT(293 

^1,1 

O.DO 

P0TENT(294-295) 

$2 1 through S 2 2 

2*0. DO 

P0TENT(296,297,...,550) 

Coefficients for Moon Model 


POTENT (296) 


O.DO 

P0TENT(297) 

^2 

2.07108D-4 

POTENT (298) 


-2.1D-5 

POTENT (299-310) 

through J-jg 

12*0. DO 

P0TENT(311) 

^1,1 

O.DO 

POTENT (31 2) 

^2,1 

O.DO 

P0TENT(313) 

^2,2 

2.0716D-5 

P0TENT(314) 

^3,1 

3.4D-5 

POTENT (31 5) 

^3,2 

O.DO 

POTENT (31 6) 

^3,3 

2.583D-6 

P0TENT(317-320) 

C^^^ through C^^^ 

4*0. DO 

P0TENT(321-325) 

Cg ^ through Cg ^ 

5*0. DO 





h 


Table 3.1. SYSTEM Namelist (Continued) 

VARIABLE DIM TYPE DESCRIPTION UNITS DEFAULT 



P0TENT(326-331) 

Ce.I S,6 

6*0. DO 

P0TENT(332-338) 

Cj ^ through Cy y 

7*0. DO 

P0TENT( 339-346) 

Cg’, through Cg’g 

8*0. DO 

POTENT (347- 355) 

Cg^l through Cg^g 

9*0. DO 

P0TENT(356-365) 

^lOJ C-jQjQ 

10*0. DO 

P0TENT(366-376) 

^11,1 C-,^ 

11*0. DO 

P0TENT( 377-388) 

C-J2 1 through C-|2 -|2 

12*0. DO 

POTENT (389-401) 

Clsil *'”'““9'' ‘=13,’l3 

13*0. DO 

P0TENT(402-415) 

'-14,1 through 

14*0. DO 

POTENT (41 6-430) 

Cisil through C,g[,g 

15*0. DO 

POTENT (431) 

^1 ,1 

O.DO ■ 

P0TENT(432-433) 

S2 1 through S2 2 

2*0. DO 

POTENT (434-436 

S^ 1 through S^ 3 

3*0. DO 

POTENT (437-440) 

S^^^ through S^^^ 

4*0, DO 

P0TENT(441-445) 

Sg 1 through Sg g 

5*0. DO 

P0TENT(446-451) 

Sg[, through Sg[g 

6*0. DO 

POTENT( 452-458) 

Sy 1 through Sy y 

7*0. DO 

POTENT (459- 466) 

Sg’, through Sg’g 

8*0. DO 

POTENT( 467-475) 

Sg^l through Sg^g 

9*0. DO 

P0TENT(476-485) 

Sio.l through sjg_,g 

10*0.00 

P0TENT(486-496) 

Sll.l through 

11*0. DO 

POTENT (497- 508) 

Si 2 1 through S-j 2 12 

12*0. DO 





Table 

3.1. SYSTEM Namelist (Concluded) 



VARIABLE 

DIM TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

POTENT( 509-521) 
POTENT (522-535) 
POTENT( 536-550) 

TEDELT 

VEDELT 

1 DP 

1 DP 

^13,1 ^13,13 

514.1 5i4.14 ■ 

515.1 

TORS ephemeris storage time delta 
Shuttle ephemeris storage time delta 

hrs 

hrs 

13*0. DO 
14*0. DO 
15*0. DO 

.045D0 

.015D0 


r\> 


3.2 MISVEH NAMELIST 

The namelist MISVEH allows the user to input all data pertaining to 
vehicle state and mission characteristics. 

3.2.1 Shuttle Vehicle State 

The input state of the Shuttle vehicle is input through the MISVEH name- 
list, and the user has the option to choose the reference frame for the vehicle 
input state. Inputs are the following: 

(a) Input state (INVEC) 

(b) Time tag of input state (HMS) 

(c) Description keys of input state (KEYS) 

(d) Indicator for transformation of input state to M50 coordinates (ITRF) 

3.2.2 TORS Vehicle State 

For OD problems which include relayed tracking data, the user may input 
the number of TORS vehicles involved and their associated input state vectors. 
Unlike the Shuttle vehicle state, TORS state vectors must be input in M50 
coordinates only. Inputs are the following: 

(a) . Number of TDRS vehicles involved (NTDRS) 

(b) Input TDRS state vectors (TDREST, TDRWST, TDRSPR) 

3.2.3 Mission History Tables 

The MISVEH namelist provides the chronological history tables for vehicle 
attitude, weight and door configuration, and nominal vent data. Inputs are 
the following: 


(a) Vehicle attitude history 

(i) Indicator to get attitude data (IBATT) 

(ii) Attitude history table (BTAB) 

(iii) Attitude rate table (ETAB) 

(b) Vehicle weight and configuration history (WTAB) 

(c) Vehicle nominal vent table (VENTIN) 

3.2.4 Vehicle Drag Parameters 

The user can set the constants for vehicle mass and cross-sectional area 
used in the vehicle drag computations. In addition, user may specify constant 
area drag by input of a nonzero drag coefficient Cg. 

Inputs are the following: 

(a) Vehicle mass (WT) 

(b) Reference cross-sectional vehicle area (CSA) 

(c) Drag Coefficient Cg (CDRAG) 

3.2.5 Mission Time Parameters 

Inputs of various time parameters associated with the vehicle mission that 
are provided in the MISVEH namelist are: 

(a) Epoch time (MONTH, DAYS, YEAR) 

( b ) Vehicle launch date -. ( GET ) 

(c) E.T. - U.T. offset calculation constants (CDETUT) 
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Table 3.2. MISVEH Namelist 


VARIABLE 



DIM 

TYPE 

DESCRIPTION 

■ UNITS 

DEFAULi’ 

BTAB 



(12,20) 

DP 

Attitude History Table 



BTAB(1,J) 

(J=l,2,.. 

ro 

o 



Entry times (in chronological order) 

hrs from 
epoch 

20*0. DO 

BTAB(2,J) 

(J=l,2,.. 

.,20) 



Vehicle attitude mode 
=1, inertial hold (IH) 

=2, inertial barbecue (BBQ) 

=3, local vertical-local horizontal (LVLH) 
=4, solar inertial (SI) 


20*0 

BTAB(I.J) 

(1=3,4,.. 

J=l,2,.. 

..11; 

.,20) 



Attitude matrix (loaded by columns) 
BTAB (2,0) = 1 =*>B0DY to M50 
BTflB(2,J) = 2=i>B0DY to M50 
BTAB (2,0) = 3=^ BODY to LVLH 
BTAB(2,0) = 4=^ BODY to SI 

■ 

180*0. DO 

BTAB(12,J) 

(J=l,2,. 

..,20) 



Pointer into ETAB for rotation information 
=N, go to ETAB(I,N) (1=1 ,2,3,4) for rotation 
axis vector and rotation rate 


20*0 

CDETUT 



2 

DP 

Coefficients in the ET-UT offset calculations 



CDETUT(l) 






sec 

21.251711D0 

CDETUT (2) 






sec/day 

.002592D0 

CDRAG 



1 

DP 

Drag coefficient C^ 

=0, variable drag model on 

if^O, value is C^ for constant drag 


O.DO 

CSA 



1 

DP 

Reference vehicle cross-sectional area 





Table 3.2. MISVEH Namelist (Continued) 


- 

VARIABLE 

DIM 

TYPE 


DESCRIPTION 

UNITS 

DEFAULT 


DAYS 


1 

R 

Day of epoch 

calendar 

day 

0 

1 

1 ETAB 


(4,5) 

DP 

Attitude 

Rate Table 



S 

i 

ETAB(I,J) 

(I=1,2,3;J=1,2,.. 

.,5) 


Axis of rotation in BODY coordinates 

ER 

3*0. DO 

f; 

J 

'i 

1 

ETAB(4,J) 

(J=l,2,...,5) 



Angular rate of rotation 

rad/sec 

O.DO 


GET 


3 

DP 

Vehicle launch date 




1 

GET(l) 




hours 


hr 

O.DO 

'!■ 

\ CO 

GET(2) 




minutes 

. from epoch 

min 

O.DO 



GET(3) 




seconds 


sec 

O.DO 


* HMS 


3 

DP 

Time tag of vehicle input state vector 





HMS(l) 




hours 


hr 

O.DO 



HMS(2) 




minutes 

. from epoch 

min 

O.DO 



HMS(3) 




seconds 

J 


sec 

O.DO 



IBATT 


1 

I 

Attitude data flag 

=0, no data 

=1, get attitude data 


0 



INVEC 


6 

DP 

Shuttle input state vector 

ER. 

ER/hr 

6*0. DO 





VARIABLE 


I TOP 


KEYS 

I KEYS(l) 


CO 

-4 


KEYS (2) 


KEYS (3) 


KEYS(4) 


Table 3.2. MISVEH Namelist (Continued) 
DIM TYPE DESCRIPTION 


DEFAULT 


1 I Indicator flag for transformation of input ER, 0 

to M50 coordinates ER/hr 

=0, no 
=1 , yes 

8 I Input state description keys 

Input state coordinates 1 

=1, x,y,z,x,y,z (CARTESIAN) 

=2, v,X,\jj,^,X,r (POLAR) 

(r is replaceable by h: See KEYS(7)) 

=3, a,e,i,fi,o),v (ORBITAL ELTS.) 

Input reference frame i 

=1, M50 
=2, M50 

=3, Selenographic (See KEYS(8)) 

=4, Skylab (not used) 

=5, MNBY 


Units for distance parameters 1 

- =1. ER 
=2, feet 
=3, yards 
=4, kilometers 
=5, nautical miles 
=6, statute miles 


Units for distance part of velocity parameters 

~1 j ER 

=2, feet 

=3, yards 

=4, kilometers 

=5, nautical miles 

=6, statute miles 


1 




Table 3.2. MISVEK Namelist (Continued) 
VARIABLE DIM TYPE DESCRIPTION 


KEYS (5) 

KEYS(6) 
KEYS (7) 


KEYS(8) 

MONTH 1 R 

NTDRS ] I 

TDREST 7 DP 

TDREST(l) 

TDREST(I) (I = 2,3,. ...7) 


Units for time part of velocity parameters 
=1 , hours 
=2, minutes 
=3, seconds 

Units for angle parameters 
=1 , radians 
=2, degrees 

Polar coordinate flag (set only in case 
KEYS(l) = 2) 

=-l, h is assumed 
=1 , r is assumed 

Selenographic reference frame flag (set only 
in case KEYS(2) = 3) 

=0, non- rotating frame 
=1 , rotating frame 

Month of epoch 


Number of TORS vehicles (NTDRS £ 3) 

TDRS-EAST state vector and time tag 
Time tag of TDRS-EAST state vector 


TDRS-EAST state vector (position and velocity) 


DP 



UNITS DEFAULT 


1 


1 


calendar 0 
month 

0 

7*0. DO 


hrs from 
epoch 

ER, 

ER/hr 


TDRSPR 


7 


TDRS-SPARE state vector and time tag 


7*0. DO 


Table 3.2. MISVEH Namelist (Continued) 


VARIABLE DIM TYPE 

TDRSPR(l) 

TDRSPR(I) (I = 2,3,. ..,7) 

TDRWST 7 DP 

TDRWST(l) 

TDRWST(I) (12,3,. ..,7) 

VENTIN (3,20) DP 

VENTIN(1,J) (J = 1,2,. ..,10) 

VENTIN(2,J) (J = 1,2,. ..,10) 

VENTIN(3,J) (J = 1,2,. ..,10) 

VENTIN(1,J) (J = 11, 12,. .,,20) 
VENTIN(2,J) (J - 11, 12,.. .,20) 
VENTIN(3,J) (J = 11,12,.. ..,20) 


DESCRIPTION 

Time tag of TDRS-SPARE state vector 

TDRS-SPARE state vector (position and velocity) 

TDRS-WEST state vector and time tag 
Time tag of TDRS-WEST state vector 

TDRS-WEST state vector (position and velocity) 

Input nominal vent table 

Vent start time (in chronological order) 

Vent stop time 

Pointer into the right half of the VENTIN 
array for vent thrust information 
=N, go to VENTIN(I,N+10) (1=1 ,2,3) for vent 
thrust values along each of the vehicle 
body axes. 

Thrust along Xgy-axis of indicated vent 
Thrust along Ygy-axis of indicated vent 
Thrust along Zgy-axis of indicated vent 


UNITS 


hrs from 
epoch 

ER, 

ER/hr 


hrs from 
epoch 


ER, 

ER/hr 


hrs from 
epoch 

hrs from 
epoch 


lbs 

lbs 

lbs 


DEFAULT 


7*0. DO 

10*0. DO 
10*0. DO 
10*0. DO 

10*0. DO 
10*0. DO 
10*0. DO 


Table 3.2. MISVEH Namelist (Concluded) 


VARIABLE 


DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

WT 



1 

DP 

Vehicle mass 

Ibs-m 

O.DO 

WTAB 



(3,5) 

DP 

Weight and Configuration History Table 



WTABd ,J) 

(0-1,2,.. 

.,5) 



Entry times (in chronological order) 

hrs from 
epoch 

5*0. DO 

WTAB(2,J) 

(0=1,2,.. 

.,5) 



Vehicle weight 

Ibs-m 

5*0. DO 

WTAB (3,0) 

(0=1,2,.. 

.,5) 



Vehicle configuration 
=1 , doors shut 
=2, doors open 


5*0. DO 

YEAR 



T 

I 

R 

Year of epoch 

calendar 

0.00 


3.3 VADRAG NAMELIST 


Whenever CDRAG = 0 in the MISVEH namelist inputs, the variable area drag 
model is called to compute the vehicle on-orbit drag coefficient C^. The 
VADRAG namelist contains the drag model parameters used to compute Cq. 

3.3.1 Drag Model Parameters 

Each parameter in the drag model is assigned two values: the first value 

is used when drag is computed for the vehicle doors-open configuration, and 
the second value is used for vehicle doors-closed configuration. Inputs to 
the drag model are: 

(a) Frontal area constant (CDF) 

(b) Side area constant (CDS) 

( c ) Normal (viewed from above) area increase constant ( CDN ) 

(d) Exponential shaping factor (CRN) 

(e) Model constant (CDA) 
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Table 3.3. VADRAG Namelist 


VARIABLE 


DIM TYPE 


DESCRIPTION 


UNITS DEFAULT 


CDA 

CDA(l) 

CDA(2) 

CDF 

CDF(l) 

CDF(2) 

CDN 

CDN(l) 

CDN(2) 

CDS 

CDS(l) 

CDS(2) 

CPN 

CPN(l) 

CPN(2) 


DP Variable drag model constant 

Doors open 
Doors closed 

DP Frontal area drag constant 

Doors open 
Doors closed 

DP Normal (viewed from above) area increase constant 

Doors open 

Doors closed 

DP Side area constant 

Doors open 
Doors closed 

DP Exponential shaping constant 

Doors open 
Doors closed 


.77015D0 

.64592D0 

.85532D0 

.7859D0 

2.74618D0 

2.41256D0 

1.84306D0 

1.92837D0 

1.31376D0 

1.31586D0 


3.4 PRBCTL NAMELIST 


The PRBCTL namelist allows the user to select the desired job mode for the 
Bench Program, and to choose from the output options available for residual 
print and display. In addition the user may select via PRBCTL the optional 
acceleration models available in the Bench Program. 

3.4.1 Problem Description 

The user has the option of providing a brief description of the current 
OD case, which description will appear in the printout summarizing the case 
input. This description may use up to 72 alphanumeric characters, including 
blanks. Input is accomplished via the PRODES array, which contains twelve 
elements, each six characters in length. 

3.4.2 DC Job Mode 

The DC job mode is input via the PRBCTL variable IDCM. User may select 

(a) S uperbatch (SB) mode (IDCM = 1) 

(b) Batch-to-^batch (BB) mode (IDCM = 2) 

(c) DC-Edit (DCE) mode (IDCM = 3) 

(d) ■ Trajectory only (TONLY) mode (IDCM = 4) 

Additional controlling input for the OD process is provided in the namelists 
SBCTL (for SB mode), BTBCTL (for BB and DCE modes), and TRAP (for TONLY mode). 

3.4.3 DC Output Options 

User option for DC output available in the PRBCTL namelist are the 
following: 

(a) Residual print option (KPR4) 

(b) DC display reference system (KREF) 
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Note that the residual print flag KPR4 serves as an indicator flag for the 
Batch Residual Plotter. See the BRESPL namelist for plotter control input. 


3.4.4 Acceleration Model Flags 

Inclusion of the drag model, solar radiation model, or vent model in the 
computations of vehicle acceleration can be accomplished via the model indi- 
cator flags provided in PRBCTL : 

(a) Drag model indicator flag (JDRAG) 

(b) Solar radiation model indicator fla g (J SOLAR) 

(c) Vent model indicator flag (JVENT) 

These models will not be considered unless the user specifies their inclusion. 


Table 3.4 


PRBCTL Namelist 


VARIABLE 


DIM TYPE 


DESCRIPTION 


UNITS 


DEFAULT 


IDCM 

JDRAG 

JSOLAR 

JVENT 

KPR4 

KREF 


1 I Job mode flag 

=1 , superbatch (SB) 

=2, batch- to-batch (BB) 

=3, DC edit (DCE) 

=4, trajectory only (TONLY) 

1 I Drag model indicator flag 

=0, model off 
=1 , model on 

1 I Solar radiation model indicator flag 
-0, model off 
=1 , model on 

1 I Vent model indicator flag 

=0, model off 
=1 , model on 

1, I Residual print option flag 

=0, print option off 
=1, print last iteration only 
=2, print first and last iteration only 
=3, print all iterations 

1 I Display reference system flag 

=1, M50 

=2, True of date 
=3, MNBY 

12 H Case description 


-1 


0 


0 


0 


3 


1 


PRODES 


12*6HJ5 


3.5 SBCTL NAMELIST 


The SBCTL namelist inputs provide the user control of the OD process in 
the superbatch operating mode. Input is required in SBCTL only in the case 
IDCM = 1 (SB mode) in PRBCTL . 

3.5.1 Batch Processing Control 

Control of data batches to processed is provided through SBCTL by 
allowing the user to specify the data batch span to be considered in the 
current superbatch OD case. The user identifies the first batch and last 
batch to be processed by their batch numbers on the Data Batch tape; 

(a) Batch number of first batch to be processed (NFBAT) 

(b) Batch number of last batch to be processed (NLBAT) 

The user is provided further control over the batch processing through the 
manual edit and batch exclusion options in the MANEDT namelist. Note: Data 

type exclusion capabilities are available in the SB mode. Parameters governing 
data type exclusion can also be found in MANEDT . 

3.5.2 Superbatch Logic Options 

Available to the user are four major options in the main superbatch logic 
of the Bench Program: 

(a) Option to anchor vehicle state at last observation time . If this 
option is desired, the user sets the logical element LANCHL = .TRUE.. 

(b) A priori covariance option . If a priori covariance is to be supplied 
in the current OD case, the user sets the logical element 

LAPCOV = .TRUE.. Covariances are input via the COVARS namelist. 


(c) Force option . If one additional iteration is desired following 
normal termination of the DC process, the user sets the logical 
element LFORCE = .TRUE.. 

(d) Option to propagate updated vehicle state after DC . If this option 
is desired, the user sets the logical element LTPROP = .TRUE.. 
Parameters required for this trajectory propagation are input via 
the TRAJ namelist. 

3.5.3 Solve-for Parameters 

In the superbatch mode, a maximum of nine solve-for parameters may be 
included in the OD solution. The user may indicate the number and type of 
solve-for parameters with the following input variables: 

(a) Solve-for relay bias flags (ISLVBI) 

(b) Solve-for drag flag (ISLVDR) 

(c) Solve-for vents flag (NVNTSL) 

A priori values for the solve-for parameters are input via the SOLFOR namelist. 

3.5.4 DC Parameters 

User may input through the SBCTL namelist all parameters associated with 
the iterated DC process. Inputs are the following: 

(a) Maximum DC iteration (MAXIT) 

(b) Maximum divergent DC iterations (MAXDV) 

(c) Minimum elevation criteria (MINEL, RMNEL, MINCA) 

(d) Convergence criteria (RCONV, VCONV; BCGNV, DCONV, VTCONV) 


3.5.5 Input State Modification 


The user may specify through the variable DA a modification to the semi- 
major axis of the vehicle input state vector. The value given to DA will be 
used to calculate new a priori position and velocity for the current OD 






Table 3.5. SBCTL Namelist 


: VARIABLE 


DIM 

TYPE 

DESCRIPTION UNITS 

DEFAULT 

i BCONV 


1 

DP 

Relay bias convergence criterion MHz 

TED 

; DA 

^ I 


1 

DP 

Delta to be added to the semi -major axis of 
input state 

0.D0 

j DCONV 


1 

DP 

Drag convergence criterion 

l.OD-2 

j ISLVBI 


9 

H 

Solve-for relay bias flags 

Flag values may be chosen in any order from the 
following list: 

9*6Ht!5 





'ESTEST', TDRS EAST 2- or 3-way relay configuration 
'WSTWST', TDRS WEST 2- or 3-way relay configuration 






'SPRSPR' , TDRS SPARE 2- or 3-way relay configuration 






'ESTWST', TDRS EAST-WEST hybrid relay configuration 


1 




'ESTSPR', TDRS East-SPARE hybrid relay configuration 
■WSTESr, TDRS WEST-EAST hybrid relay configuration 
'WSTSPR' , sDRS WEST-SPARE hybrid relay configuration 


'■ 1 




'SPREST', TDRS SPARE-EAST hybrid relay configuration 
'SPRWST', TDRS SPARE-WEST hybrid relay configuration 


} ISLVDR 


1 

I 

Solve-for drag flag 
=0, no solve-for drag 
=1, solve for drag 

0 

LANCHL 


1 

L 

Logical element 

=.TRUE., anchor vehicle state at last observation 

.FALSE. 





time 


LAPCOV 


1 

L 

Logical element 

=.TRUE., a priori covariance supplied 

.FALSE. 

L FORCE 


1 

L 

Logical element 
=.TRUE., force option on 

.FALSE. 



Table 3.5. SBCTL Namelist (Concluded) 

VARIABLE 

DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

LTPROP 

1 

L 

Logical element 

=.TRUE., trajectory propagation after DC is desired 

.FALSE. 

MAXDV 

1 

I 

Maximum divergent iterations in DC 


1 

MAXIT 

1 

I 

Maximum iterations in DC 


1 

MINCA 

1 

DP 

Minimum altitude of closest approach for TORS/ 
Shuttle RF link 

ER 

2.7437444445 

MINEL 

1 

DP 

Minimum direct elevation 

rad 

.87266463D-1 

NFBAT 

1 

I 

Number of first batch to process 


-1 

cn 

o NLBAT 

1 

I 

Number of last batch to process 


-1 

NVNTSL 

1 

I 

Number of solve-for vents 

=0, no solve-for vents 

>0, solve for N vents (where NVNTSL = N) 


0 

RCONV 

1 

DP 

Position convergence criterion 

ER 

l.OD-8 

RMNEL 

1 

DP 

Minimum relay elevation angle 

rad 

.87266463D-1 

VCONV 

1 

DP 

Velocity convergence criterion 

ER/hr 

l.OD-5 

VTCONV 

1 

DP 

Solve-for vent convergence criterion 

lbs 

l.OD-4 


3.6 BTBCTL NAMELIST 


The BTBCTL namelist inputs provide the user contr '' of the OD process in 
i-'.e batch-to-batch operating modes. Input is required in BTBCTL only in the 
cr;-. e IDCM = 2 (BB mode) or IDCM = 3 (DCE mode) in PRBCTL . 

3.6.1 Batch Processing Control 

Control of data batches to be processed is provided by allowing the user 
to specify the data batch span to be considered in the current batch-to-batch 
OD case. User identifies the first batch and last batch to be processed by 
their batch numbers on the Data Batch Tape: 

(a) Batch number of first batch to be processed (NFBAT) 

(1^) Batch number of last batch to be processed (NLBAT) 

The user is provided further control over the batch processing through the 
manual edit and batch exclusion options in the MANEDT namelist. 

3.6.2 Batch-to-Batch Logic Options 

Available to the user are three major options in the main batch-to-batch 
logic of the Bench Program: 

(a) A priori covariarice option . If a priori covariance is to be supplied 
in the current OD case, the user sets the logical element 

LAPCOV = .TRUE.. Covariances are input via the namelist COVARS. 

(b) Automatic edit option . If the Batch Edit option is desired in the 

BB mode (IDCM = 2), the user sets the logical element LEDIT = .TRUE.. 
Parameters for the automatic editor are available in the ATOEDT 
namelist. Note: If the DCE mode has been selected (IDCM = 3), it 

is mandatory that the user set LEDIT = .TRUE.. 


(c) Force option . The force option, enabled by the logical element LFORCE, 
is available in the batch-to-batch operating modes. If a data batch 
is to be processed in one of these modes, and an additional DC 
iteration is to be forced following normal termination of the DC 
process in that mode, the user sets LFORCE = .TRUE.. If the force 
option is called in either the BB mode with edit option, or in the 
DCE mode, no editing is performed subsequent to the forced iteration. 

Note : Care should be taken with the use of this option in the 

batch-to-batch modes. The force option will apply to all data batches 
processed in the batch-to-batch modes. Therefore, if only a partic- 
ular data batch is to be forced, that batch should be processed as 
a single batch case. 

3.6.3 S olve-for Parameters 

Solve-for capability in the BB or DCE mode is limited to a single solve- 
for vent. Indication of a solve-for vent is accomplished with the solve-for 
vent flag NVNTSL. 

Note: No a priori vent information is input in the batch-to-batch modes. 

A priori vent force levels are assumed to be zero, and the vent time span is 
assumed to be that of the data batch. 

3.6.4 DC Parameters 

User may input via the BTBCTL namelist al^ parameters associated with the 
DC iterative process. Inputs are the following: 

(a) Maximum DC iterations (MAXIT) 

(b) M aximum divergent DC iterations (MAXDV) 


(c) Minimum direct elevation angle (MINEL) 

(d) Minimum relay elevation angle (RMNEL) 

(e) Minimum altitude for TDRS/Shuttle RF link (MINCA) 

(f) Convergence criteria (RCONV, VCONV: VTCONV) 


Table 3.6. BTBCTL Namelist 


VARIABLE 

LAPCOV 

LEDIT 

LFORCE 

MAXDV 
MAXIT 
cn MINCA 

-p» 

MINEL 

NFBAT 

NLBAT 

NVNTSL 

RCONV 

RMNEL 

VCONV 


DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

1 

L 

Logical element 

=.TRUE., a priori covariance supplied 


.FALSE. 

1 

L 

Logical element 

=.TRUE., automatic edit desired 


.FALSE. 

1 

L 

Logical element 
=.TRUE., force option on 


.FALSE. 

1 

I 

Maximum divergent iterations in DC 


1 

1 

I 

Maximum iterations in DC 


1 

1 

DP 

Minimum altitude of closest approach for 
TDRS/Shuttle RF link 

ER 

2.7437444445D 

1 

DP 

Minimum direct elevation 

rad 

.872664630-1 

1 

I 

Number of first batch to process 


-1 

1 

I 

Number of last batch to process 


-1 

■; 

I 

Number of solve-for vents 
=0, no solve-for vents 
=1, solve for single vent 


0 

1 

DP 

Position convergence criterion 

ER 

l.OD-8 

1 

DP 

Minimum relay elevation angle 

rad 

.87266463D-1 

1 

DP 

Velocity convergence criterion 

ER/hr 

l.OD-5 

1 

DP 

Solve-for vent convergence criterion 

lbs 

l.OD-4 


VTCONV 


3.7 TRAJ NAMELIST 


The TRAJ namelist variables are associated with the control of trajectory 
and trajectory-related output. Input in the TRAJ namelist is required when 
a propagation of vehicle state is called for in the TONLY mode (IDCM = 4), 
or in the SB mode with the LTPROP option (propagation of updated state). 

3.7.1 Trajectory Parameters 

Trajectory-related inputs in the TRAJ namelist are the following: 

(a) Trajectory start time (TSTART) 

(b) Trajectory stop time (TSTOP) 

(c) Trajectory printout time delta (DELTAT) 







DELTA! 

1 

DP 

Trajectory printout time delta 

hrs 

1.0 

TSTART 

1 

DP 

Propagation start time 

hrs from 
epoch 

O.DO 

TSTOP 

1 

DP 

Propagation stop time 

hrs from 
epoch 

O.DO 
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3.8 SOLFOR NAMELIST 


The SOLFOR namelist allows the user to input a priori values for the 
solve-for parameters selected for the current OD case. Input is required 
for the SOLFOR namelist only in the case that solve-for parameters are in- 
cluded in a superbatch OD solution. Note: A batch-to-batch OD solution in 

which a vent is included does not require SOLFOR input. In that setting, 
a priori vent force levels are assumed to be zero, and vent start and stop 
times are identified with the times of the first and last valid observations, 
respectively, of each data batch. 

3.8.1 A Priori Solve-for Parameters 

A priori solve-for inputs are the following: 

( a ) A priori solve-for vent inputs 

(i) A priori vent force levels (VTFORC) 

(ii) Solve-for vent start and stop times (VTSTRT, VTSTOP) 

(b) A priori solve-for drag multiplier (DKFACT) 

(c) A priori solve-for relay bias values (BIAS) 







Table 3.8. SOLFOR Namelist 




VARIABLE 

DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

BIAS 

9 

DP 

A priori relay bias values (in order of bias 
flags in ISLVBI array) 

MHz 

9*0. DO 

DKFACT 

1 

DP 

A priori drag multiplier 


1.0 

VTFORC 

(3,3) 

DP 

A priori solve-for vent forces 




VTF0RC(1,J)(J=1,2,3) 



Thrust along Xgy-axis of J-th vent 

lbs 

3*0. DO 


VTF0RC(2,B) (J=l ,2,3) 



Thrust along Ygy-axis of J-th vent 

lbs 

3*0. DO 


VTF0RC(3,J) (J=l,2,3) 



Thrust along Zgy-axis of J-th vent 

lbs 

3*0. DO 

07 VTSTRT 
00 

3 

DP 

Solve-for vent start times 

(in order of vents in VTFORC array) 

hrs from 
epoch 

3*0. DO 

VTSTOP 

3 

DP 

Solve-for vent stop times 

(in order of vents in VTFORC array) 

hrs from 
epoch 

3*0. DO 


3.9 COVARS NAMELIST 


User input of a priori covariances and associated parameters is accom- 
plished via the COVARS namelist. Input is required only in the case that the 
a priori covariance option is selected by the user (LAPCOV = .TRUE, in SBCTL 
or BTBCTL ). 

3.9.1 State Covariance Reference Frames 

The user may select the reference frames both for input and display of the 
vehicle state covariance: 

(a) State covariance input reference frame (ICOV 

(b) State covariance display reference frame (IDIS) 

3.9.2 Covariance Matrix Elements 

The user may specify the elements of the a priori covariance matrix with 
with the following inputs: 

(a) A priori state covariance (SSCOV) 

(b) A priori vent covariance (SCVNT) 

(c) A priori drag covariance (SDRAG) 

(d) A priori relay 2- or 3-v/ay relay Doppler bias covariance (SCVB23) 

(e) A priori relay hybrid Doppler bias covariance (SCVBH) 

The actual a priori covariance matrix constructed from these elements will be 
a block diagonal matrix A of maximum dimension 15x15. The blocks on the diagonal 
are input by the user from the above list, and the order of the diagonal follows 
the solution vector order: 


sscov 



Note: In the case of multiple vents specified in the solution vector, SCVNT 

is used as the a priori covariance for each vent. 

3.9.3 KGAMMA Downweighting 

The user may take the option to downweight the constructed a priori 
covariance matrix (KGAMMA downweight option). Parameters associated with this 
option are: 

(a) KGAMMA covariance multiplier (FKGAMMA) 

(b) Number of times KGAMMA PBI is depressed (NP9I) 


Table 3.9. COVARS Namelist 


VARIABLE 

DIM 

TYPE 

DESCRIPTION UNITS 

DEFAULT 

FKGAMA 

1 

R 

K6AMMA covariance multiplier 

2.0 

ICOV 

1 

I 

Input reference frame for state covariance 

=1, M50 coordinates 

=2, UVW coordinates 

=3, classical orbital elements 

1 

IDIS 

1 

I 

Display reference frame for state covariance 
=0, M50 coordinates only 
=1, M50 coordinates, UVW coordinates, and 
classical orbital elements 

0 

NPBI 

1 

I 

Number of times Kgamma PBI is depressed 

0 

SCVB23 

1 

DP 

A priori relay 2- or 3-way Doppler bias covariance 

O.DO 

SCVBH 

1 

DP 

A priori relay hybrid Doppler bias covariance 


SCVNT 

6 

DP 

A priori vent covariance matrix 
(lower triangular portion, by rows) 

6*0. DO 

SDRAG 

1 

DP 

A priori drag covariance 

O.DO 

SSCOV 

21 

DP 

A priori state covariance matrix 
(lower triangular portion, by rows) 

21*0. DO 


3.10 ATOEDT NAMELIST 


The ATOEDT namelist allows the user to set the parameters governing the 
operation of the automatic editor. The editor is called in the DCE mode 
and in the BB mode with edit option on. 

The edit function sequence is as follows: gross pre-edit, D1 (divided 

first differences) edit, D2 (divided second differences) edit, re-edit, and 
edit state testing. The user may set the controlling parameters for each 
of these edit functions. In addition, the user specifies the number of times 
the editor is called by the OD executive in the Bench Program. 

3.10.1 Edit Loop Limit 

An edit loop consists of an edit pass of the current residuals of a data 
batch followed by a DC on the edited data batch. The number of edit loops 
called for is set by the user via the MXELOP variable. 

3.10.2 Gross Pre-edit Parameters 

Parameters for the gross pre-edit are the following: 

(a) Data type edit limits (RMAX) 

(b) Percentage of edited points allowed (EDLIM) 

3.10.3 D1 Edit Parameters 

Parameters for the D1 edit are the following: 

(a) Minimum number of points for D1 edit (MINI) 

(b) Minimum editing criteria (RMIN) 

(c) D1 edit limit search parameters (DELTAl, DKl , SRCHl) 


3.‘;J.4 D2 Edit Parameters 


Parameters for the D2 edit are the following: 

(a) Minimum number of points for D2 edit (MIN2) 

(b) Minimum editing criteria (RMIN) 

(c) D2 edit limit search parameters (DELTA2, DK2, SPCH2) 

3.10.5 Edit State Testing Parameters 

In the edit state tests, each data type is tested for 1) the presence of 
large data jumps, 2) the presence of large groups of consecutive edited points, 
and 3) an excessive total number of edited points. Parameters for the edit 
state tests are the following: 

(a) Data discontinuity multiplier (ID) 

(b) Maximum number of consecutive edited points (IG) 

(c) Percentage of edited points allowed (EDLIM) 
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Table 3.10. ATOEDT Namelist 


VARIABLE 

DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

DELTAl 

1 

DP 

D1 edit parameter which determines the largest 
jump between successive Dlb values 


.5 

DELTA2 

1 

DP 

D2 edit parameter which determines the largest 
jump between successive D2b values 


1.5 

DKl 

1 

DP 

D1 edit parameter which determines the largest jump 
of Dlb values relative to the reference Dlb. 

4.0 

DK2 

1 

DP 

D2 edit parameter which determines the largest 
jump of D2b values relative to the reference 
D2b. 


12.0 

EDLIM 

1 

R 

Percentage of edited points allowed 
(pre-edit, edit state tests) 


.8 

ID 

1 

I 

Data discontinuity multiplier 
(edit state tests) 


100 

IG 

1 

I 

Maximum number of consecutive edited points 
(edit state tests) 


100 

MINI 

1 

I 

Minimum number of points required for D1 edit 


3 

MIN2 

1 

I 

Minimum number of points required for D2 edit 


5 

MXELOP 

1 

I 

Maximum number of edit loops 


3 

RMAX 

3 

DP 

Edit limits for gross pre-edit 



RMAX(l) 



Angles edit limit 

rads 

l.ODO 

RMAX (2) 



Range edit limit 

hrs 

l.ODO 

RMAX(3) 



Doppler edit limit 

MHz 

l.ODO 


Table 3.10. AlOEDT Namelist (Concluded) 


VARIABLE 

DIM 

TYPE 

DESCRIPTION 


UNITS 

DEFAULT 

RMIN 

3 

DP 

Minimum editing criteria for D1 

and D2 edit 



RMIN(l) 



Angles edit minimum 


rads 

O.DO 

RMIN(2) 



Range edit minimum 


hrs 

O.DO 

RMIN(3) 



Doppler edit minimum 


MHz 

O.DO 

SRCHl 

1 

R 

Level-setting parameter for Dlb 
value 

reference 


.7 

SRCH2 

1 

R 

Level -setting parameter for D2b 
value 

reference 


.7 


3.11 MANEDT NAMELIST 


The MANEDT namelist furnishes to the user all manual editing and data 
exclusion options available in the Bench Program. Some options are available 
in superbatch mode only, and will be indicated by (SB only). 

3.11.1 Manual Data Editing Instructions 

The user may input all data necessary to perform manual editing via the 

variable array MEDIT. The MEDIT inputs allow the user to edit batch data by 

frames. 

3.11.2 Data type exclusion Instructions (SB only) 

In the SB mode, the user has the option to exclude specific data types 
from all batches to be processed. Each permissible data type exclusion has 
its own controlling logical element. If the user desires a certain data type 
to be excluded, he sets the appropriate logical variable to .TRUE. For 
example, exclusion of all C-band azimuth data will require the logical element 

LEXCCA = .TRUE.. See the MANEDT namelist table for a complete listing. 

3.11.3 Batch Exclusion Instructions 

Two methods are available for exclusion of data batches: 

, (a) Batch exclusion by batch number. The variable array I EXBAT can be 
assigned the batch numbers of data batches to be excluded, up to a 
maximum of 15 batches. 

(b) Batch exclusion by station type (SB only). This batch exclusion 
capability, available in superbatch mode only, is provided to the 
user for ease in setting up a superbatch OD case in the Bench 
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Program. Specifically, the user may exclude all C-band, all S-band, 
or all relay batches from those to be processed. Each type of batch 
exclusion has its own controlling logical element. Setting the 
appropriate logical element to .TRUE, will cause the batches of the 
selected type to be excluded. Logical elements are the following: 

(i) LEXCLC - Exclude all C-band batches. 

(ii) LEXCLS - Exclude all S-band jatches. 

(iii) LEXCLR - Exclude all relay batches. 


Table 3.11. MANEDT Namelist 


VARIABLE 

DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

I EXBAT 

15 

I 

Numbers of batches to be excluded 


15*0 

LEXCCA 

1 

L 

Logical element 

=.TRUE., exclude azimuth angle in all C-faand 
batches 


.FALSE. 

LEXCCE 

1 

L 

Logical element 

=,TRUE., exclude elevation angle in all C-band 
batches 


.FALSE. 

LEXCCR 

1 

L 

Logical element 

= .TRUE., exclude, range in all C-band batches 


.FALSE. 

LEXCHD 

1 

L 

Logical element 

=.TRUE., exclude hybrid Doppler in all relay 
batches 


.FALSE. 

LEXCHR 

1 

L 

Logical element 

=.TRUE., exclude hybrid range in all relay 
batches 


.FALSE. 

LEXCLC 

1 

L 

Logical element 

=.TRUE., exclude all C-band batches 


.FALSE. 

LEXCLR 

1 

L 

Logical element 

=.TRUE., exclude all relay batches 


.FALSE. 

LEXCLS 

1 

L 

Logical element 

=.TRUE., exclude all S-band batches 


.FALSE. 

LEXCRD 

1 

L 

Logical element 

=.TRUE., exclude 2- or 3-way Doppler in all 
relay batches 


.FALSE. 











Table 3.11. MANEDT Namelist (Concluded) 

VARIABLE DIM TYPE DESCRIPTION UNITS DEFAULT 


LEXCRR 

LEXCSX 

LEXCSY 

LEXCSR 

LEXCS2 

LEXCS3 

MEDIT 

MEDIT(1,J) (0=1,2 
MEDIT(2,J) (0=1,2 

MEDIT(3,0) (0=1,2 
MEDIT(4,0) (0=1,2 


1 L 

1 L 

1 L 

1 L 

1 L 

1 L 

(4,100) I 

..., 100 ) 

..., 100 ) 

..., 100 ) 

..., 100 ) 


Logical element 

=.TRUE., exclude 2- or 3-way range in all 
relay batches 

Logical element 

=.TRUE., exclude X angle in all S-band batches 
Logical element 

=.TRUE., exclude Y angle in all S-band batches 
Logical element 

=.TRUE., exclude range in all S-band batches 
Logical element 

=,TRUE., exclude 2-way Doppler in all S-band 
batches 

Logical element 

=.TRUE., exclude 3-way Doppler in all S-band 
Batches 

Manual data editing instructions 

Numbers of batches to be edited 

Data type(s) to be edited 

=0, all data types 

=1 , angle 1 

=2, angle 2 

=3, range 

=4, Doppler 

Beginning frame number for editing 
Ending frame number for editing 


.FALSE. 


.FALSE. 

.FALSE. 

.FALSE. 

.FALSE. 

.FALSE. 

100 *(- 1 ) 

100 *(- 1 ) 


100 *(- 1 ) 

100 *(- 1 ) 


3.12 BRESPL NAMELIST 


The BRESPL namelist allows the user to specify the operating parameters 
for the Batch Residual Plotter. The plotter is activated by the residual print 
control flag KPR4 in the PRBCTL namelist. 

3.12.1 Residual Plotter Units 

The user specifies the units for the residuals of each data type. This is 
done with the I UN ITS array. 

3.12.2 Residual Plotter Maximum y. Vflues 

The maximum scale values for each data type to be used by the plotter are 
input through the variable array RPLOT. Note; Default values for RPLOT are 
4xRMS of each data type. 





Table 3.12. BRESPL Namelist 


VARIABLE 

DIM 

TYPE 

DESCRIPTION 

UNITS 

DEFAULT 

I UN ITS 

3 

H 

Names of residual plot units 


■ ■ ■ ' — ' ■ 

lUNITS(l) 



Angle units 

Options are 'MILR.AD', 'RAD', or 'DEG' 


3HDEG 

IUNITS(2) 



Range units 

Options are TTPRS’, 'YDS', or 'FT' 


2HFT 

IUNITS(3) 



Doppler units 

Options are 'MHZ' or 'HZ' 


3HMHZ 

RPLOT 

4 

DP 

Residual plot maximum scale values 


■ 

RPLOT(l) 



Angle 1 maximum scale value 
=0, RPLOT(l) = 4*RMS 

input by 
lUNITS(l) 

0 . DO h 

r. 

j 

RPLOT (2) 



Angle 2 maximum scale value 
=0, RPL0T(2) = 4*RMS 

input by 
lUNITS(l) 

O.DO 1 

i- 

RPL0T(3) 



Range maximum scale value 
=0, RPL0T(3) = 4*RMS 

input by 
IUNITS(2) 

O.DO J 

RPL0T(4) 



Doppler maximum scale value 
=0, RPL0T(4) - 4*RMS 

input by 
IUNITS(3) 

O.DO 



1 


iricW.'. 


3.13 BIASCT NAMELIST 


The user has the option to apply observation biases to the data types in 
the batches to be processed in the current OD case. 

3.13.1 Observation Data Bias Instructions 

The data types to be biased and the batches in which the bias is to be 
applied are the information input via the variable array BIASIN. Bias may 
be applied to a data type in a single batch, to a data type in all batches from 
a particular ground station, or to a data type in all relayed batches having a 
particular forward or return link ID. 

3.13.2 Observation Data Bias Values 

The observation biases to be applied to the data types indicated in BIASIN 
are input via the variable array BIASV. 


Table 3.13. BIASCT Namelist 


VARIABLE DIM TYPE 


DESCRIPTION 


UNITS DEFAULT 


BIASIN (5,20) H,I 

I BIASIN(1,J) (J=l,2,...,20) H 


BIASIN(2,J) (J=l,2,...,20) I 

BIASIN(3,J) (J=l ,2,.. . ,20) H 
BIASIN(4,J) (J=l,2,...,20) H 
BIASIN(5,J) (J=l,2,...,20) H 


Observation data bias instructions 

Data type bias flags 
Flag options are: 

'SX' , S-band X angle 

'SY' , S-band Y angle 

’SR' , S-band range 

'S2D', S-band 2-vjay Doppler 

'S3D', S-band 3-way Doppler 

'A' , C-faand azimuth 

'E' , C-band elevation 

'CR' , C-band range 

'R2R', Relay 2- or 3-way range 

'R2D', Relay 2- or 3-way Doppler 

'RHR', Relay hybrid range 

'RHD', Relay hybrid Doppler 

Number of data batch in which J-th data type 
is to be biased, 

(If BIASIN(2,J) = 0, the J-th data type will 
be biased in all batches from the station or 
TORS vehicle indicated below.) 

Station ID 

(All data batches frcm this station will have 
the J-th type biased.) 

Forward link TORS ID. 

(All relayed batches with this forward link 
ID will have the J-th data type biased.) 

Return link TORS ID. 

(All relayed batches with this return link 
ID will have the J-th data type biased.) 


20*6Hb 


20*0 


20*0 


20*3H*5 

20*3H}j 



Table 3.13. BIASCT Namelist (Concluded) 


VARIABLE DIM TYPE DESCRIPTION UNITS DEFAULT 


BIASV 20 DP Observation bias values. rads, hrs, 20*0. DO 

(The value prescribed in BIASV(J) or MHz 

will be that bias applied to the J-th 
data type.) 


45 ^ 


4. PROGRAM RUN PROCEDURE 




■I 




This section gives the logical unit assignments fus; the OPS Bench Program, 
and provides the program deck setup for a Bench Program run. 

4.1 TAPE UNITS 

The OPS Bench Program requires the following two tapes for any program run: 

(a) Program Tape 

(b) Planetary Ephemeris Tape 

The current ephemeris tape in use is the JPL Tape PE50D. The PE50D tape is 
assigned to logical unit 13. 

There are two optional tapes: 

(a) Observation Data Batch Tape 

(b) Station Characteristics Tape 

These tapes are required input for a batch-to-batch or superbatch OD case, 
and may be omitted for a trajectory only case. The Data Batch Tape is 
assigned to logical unit 9; the Station Characteristics Tape is assigned to 
logical unit 10. 


4.2 PROGRAM DECK SETUP 


Figure 4-1 on the following page provides an illustration of the program 
deck setup for an OPS Bench Program run. 


















5. PROGRAM OUTPUT; SAMPLE CASES (TBD) 


This section is intended to contain several sample cases which illustrate 
typical OPS Bench Program runs. Each case will include the following: 

(a) Case description 

(b) List of program options employed in that case 

(c) Listing of card input for that case 

(d) Output print generated in the execution of that case 

The user will find in these sample cases all the types of printed output availabl 
from the Bench Program: printout of card input and data batches, DC printout, 

trajectory printout, and the various elements of display printout. 


